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ABSTRACT 

AN  INVESTIGATION  OF  THE  SIDE  THRUST  AVAILABLE  FROM 
PROPELLER  SHROUDING  SURFACES 

by 
SAMUEL  JAMES  GORDON 

Submitted  to  the  Department  of  Mechanical  Engineering  and  to  the  Department 
of  Naval  Architecture  and  Marine  Engineering  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Master  of  Science  in  Mechanical  Engineering 
and  the  Professional  Degree  of  Naval  Engineer.  . 

A  theory  and  analytical  method  are  developed  for  the  determination  of  the 
circumferential  circulation  distribution,  over  an  annular  airfoil  of  arbitrary 
included  angle.    The  resultant  circulation  distribution  is  then  used  to  determine 
the  forces  acting  on  the  airfoil  surface. 

The  theory  is  based  on  the  "L-Method"  developed  by  J.  Weissinger,  which 
consists  of  replacing  the  lifting  surface  by  a  single  concentrated  ring  vortex  and 
an  associated  system  of  trailing  vortices.    The  method  presented  here  is  com- 
pared to  other  theories  in  the  limiting  case  of  the  complete  symmetrical  annulus. 
It  is  shown  to  give  the  same  result  as  Weissinger  found  using  the  "L-Method.  " 
This,  in  turn,  is  shown  to  be  a  very  close  approximation  to  the  more  accurate 
theories  and  to  experimental  results. 

A  sample  calculation  is  made,  based  on  the  method  presented,  of  the  forces 
acting  on  a  semicircular  propeller  shroud,  to  determine  its  usefulness  as  a 
steering  device  in  lieu  of  the  conventional  ship's  rudder.    The  results  are  com- 
pared with  those  for  a  fully  movable  spade  type  rudder  with  area  equal  to  that 
of  the  shroud  and  with  a  high  lift  to  drag  ratio.     It  is  found  that  the  predicted  side 
thrust  from  the  shroud  is  more  than  twice  the  maximum  predicted  for  the  rudder. 
Additionally  the  shroud  is  found  to  give  a  thrust  augmentation  equal  to  approxi- 
mately 10%  of  the  propeller  thrust  while  the  rudder  absorbs  a  minimum  of  1% 
of  the  propeller  thrust  and  a  maximum  of  30%. 


Thesis  Supervisors:   Justin  E.  Kerwin 

Associate  Professor  of 
Naval  Architecture 

Hal  L.   Moses 
Assistant  Professor  of 
Mechanical  Engineering 
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Introduction: 

The  literature  contains  a  good  deal  of  material  on  the  theory  of  annular 
airfoils  and  ducted  propellers.    Authors  such  as  Morgan^  '  and  Weissinger'  ' 
have  treated  the  problem  of  a  symmetrical  ring  airfoil  at  an  angle  attack, 
alone,  as  well  as  in  the  presence  of  other  objects  such  as  propellers  and  hubs. 
Weissinger  has  also  treated  the  case  of  a  complete  annular  airfoil  with  deflected 
control  surfaces.    There  seems  to  be  no  evidence  of  any  previous  work  on  in- 
complete rings. 

Since  the  incomplete  shroud  seemed,  from  an  intuitive  standpoint,  to 
provide  a  very  effective  means  of  steering  a  ship,  I  constructed  a  small  self- 
propelled  model  and  fitted  it  with  a  semi-circular  propeller  shroud.    The 
model  is  illustrated  in  Figure  1.    By  rotating  the  shroud  about  the  propeller 
axis  the  side  thrust  developed  by  the  shroud  could  be  directed  so  as  to  cause 
the  model  to  turn  to  the  right  or  left.    The  steering  device  worked  very  well  on 
the  model  and  provided  the  stimulus  for  this  investigation. 

The  theory  developed  here,  I  feel,  provides  a  means  for  making  a  reasonably 
accurate  prediction  of  the  performance  of  either  complete  or  incomplete  shrouds 
with  a  circumferential  variation  in  circulation  distribution. 
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Conclusions: 

The  results  of  this  sample  calculation  indicate  that  partial  shrouds  or 
complete  shrouds  with  controllable  circulation  and  lift  distributions  can  be  used 
to  good  advantage  as  steering  devices  for  fluid  born  vehicles  such  as  ships  and 
ground  effect  machines. 

The  following  improvements  in  performance  over  that  obtained  with  con- 
ventional rudders  are  predicted: 

1.  Substantial  elimination  of  rudder  drag  at  all  turning  rates. 

2.  Thrust  augmentation  from  the  shroud. 

3.  Improved  steering  capability  with  little  or  no  way  on. 

4.  Significant  decrease  in  turning  diameter  with  an  equivalent  size 
turning  device. 

The  following  possible  shortcomings  in  the  application  of  shrouds  as 
steering  devices  require  consideration  and  investigation: 

1.  Decrease  in  propeller  loading  necessitated  by  shroud  induced 
cavitation. 

2.  Propeller  induced  shroud  vibration. 

3.  Possibility  of  loss  of  steering  capability  with  engines  stopped. 

It  is  hoped  that  further  theoretical  work  in  this  area  along  with  experi- 
mental testing  will  be  carried  out  and  that  the  apparent  benefits  to  be  gained 
from  the  use  of  propeller  shrouds  for  steering  will  be  employed  in  practice. 


SIDE  VIEW  OF  MODEL 


VIEW  SHOWING  SHROUD  IN  POSITION  FOR  TURNING 

TO  THE  RIGHT 


FIGURE  I 
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Procedure:    The  method  utilized  for  determining  the  forces  on  the  shroud  is  that 
developed  by  J.  Weissinger  and  known  as  the  "L- Method"  or  "Weis singer's 
Method.  "   The  following  basic  assumptions  are  made: 

(a)  The  fluid  is  incompressible  and  baratropic. 

(b)  Body  forces,  such  as  gravity,  are  neglected. 

(c)  The  flow  is  steady. 

(d)  Viscous  effects  are  confined  to  a  negligibly  thin  boundary  layer 
which  does  not  separate  on  the  shroud. 

The  method  consists  of  the  replacement  of  the  lifting  surface  by  a  single,  con- 
centrated bound  vortex  at  the  1/4  chord  point  and  an  associated  system  of 
trailing  vortices  as  shown  in  Fig.  2.     The  circulation  distribution,  along  the 
bound  vortex,  is  determined  by  requiring  that  the  velocity  normal  to  the  mean 
line  at  the  3/4  chord  point  must  be  zero.    Since  the  circulation  distribution  varies 
along  the  span  of  the  shroud,  a  system  of  trailing  vortices  will  be  shed  off  which 
will  form  a  cylindrical  sheet  of  radius  equal    to  that  of  the  concentrated  vortex 
ring.    The  boundary  condition  is  linearized  by  requiring  that  the  quotient  of  the 
radial  velocity  and  the  free  stream  velocity,  at  the  3/4  chord  point,  along  this 
cylindrical  sheet,  be  equal  to  the  slope  of  the  section  at  the  3/4  point.    That  is: 

PL:     — «r[cJ(i,.\3/4)-a(§)] 

?  Cl.fc\8/4> 

Once  the  circulation  distribution  is  known,  the  radial  force  and  axial  forces 
per  unit  length  at  any  section  can  be  found  by  the  "Kutta-Joukowski  Law"  which, 
in  this  case,  states: 

P2:     Fr(0)  =  -pVxptt,  e,  o)r(9) 

P3:      Fx(0)  =   p[qr    (1,  0,  0)  +  qrg(l,  9,  *l/4>]r(e) 


-  X 


TRAILING-  VORTEX 


(Thickness  not  shown) 
FIGURE  2 
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Shroud  length. 

Limiting  half  angle  of  shroud. 

Angular  coordinate  along  concentrated  vortex  in  1/4  chord  plane. 

Radius  of  concentrated  vortex  in  1/4  chord  plane. 

Pivotal  point  at  radius  equal  to  r^M  in  3/4  chord  plane. 

Angular  coordinate  of  pivotal  point  in  3/4  chord  plane. 

Displacement  of  mean  line  relative  to  nose  tail  line,  nondimensionalized 

on  L. 

Slope  of  mean  line  relative  to  nose  tail  line. 

x/ri/4  =  Nondimensionalized  axial  coordinate 

r/rl/4  =  Nondimensionalized  radial  coordinate 

Convergence  angles  of  nose-tail  cone 

Angle  of  attack  in  plane  of  symmetry. 
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For  a  shroud  which  is  symmetrical  about  8  =  0,  the  resultant  radial  force 

— ♦ 

R  lies  in  the  plane  of  symmetry  and  is  equal  to  the  total  lift  acting  on  the  shroud. 

— »  — » 

R  can  be  found  by  integrating  the  components  of  Fr(9)»  parallel  to  the  plane  of 

symmetry,  over  the  span  as  follows: 

0 


P4:      T(0)  =  V,  ..         Fr(0)  cos  9  d0 


Similarly,  the  total  induced  drag  D,  acting  on  the  shroud,  is  found  as 
follows: 

e 

^     Pi-V1/4J    F^9)d8 

The  problem  reduces  to  the  determination  of  the  circulation  distribution 
for  a  given  propeller  geometry  and  loading  and  a  chosen  shroud  geometry. 
The  circulation  distribution  is  represented  as  a  Fourier  cosine  series. 

oo 

A-23:     r(Q)  =    £   Ancos  (n0) 
n=0 

It  is  to  be  noted  here  that  9  is  the  actual  polar  angle  measured  from  the  plane 
of  symmetry.     For  an  incomplete  shroud  the  condition  that  the  circulation 
vanish  at  the  tips  is  satisfied  by  requiring  that  the  sum  of  all  terms  in  the  series 
vanish  rather  than  that  each  term  vanish  individually.    The  resultant  radial 
velocity  induced  by  the  shroud,  at  a  point  on  the  trailing  vortex  sheet,  is  deter- 
mined as  a  series  equation  in  terms  of  a  set  of  coefficients  and  integrals. 

~      6 
A~26:     5rs  =  r  £     £     C^I^ 
n=0  m=l 
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where  Cnm  are  coefficients  involving  the  parameters  An  and  Inm  are  geometric 
integrals. 

Since  four  terms  in  the  series  representation  of  F/q\  (Eq.  A-23)  will 
provide  a  representation  of  more  than  sufficient  accuracy  for  our  purposes,  we 
will  restrict  our  consideration  to  n  =  0,  . . . ,  3.    We  must,  therefore,  evaluate  a 
maximum  of  24  coefficients,  Cmn,  and  24  integrals,  Imn  for  each  angular  co- 
ordinate  ($)  at  which  we  desire  to  determine  qr  .    These  can  then  be  combined, 
using  Eq.  A-26,  to  form  an  equation  in  terms  of  the  coefficients,  An,  of  Eq. 
A-23.    If  we  let  n  =  0,  . . . ,  i^^,  nmax  ^  3,  we  must  evaluate  qr    at  (nmax  -  1) 
values  of  $  by  using  Eq.  A-26. 

Knowing  the  propeller  velocity  field  for  a  given  propeller  geometry  and 

loading,  and  knowing  the  shroud  geometry  we  can  evaluate  qr    at  the  same 

(nmax  -  1)  values  of  $  by  Eq.  P-l.     In  the  case  of  an  incomplete  shroud  we  have 

an  additional  boundary  condition  to  satisfy  and  that  is  that  IV  @\  =  0.    This  gives 

us  an  additional  equation  resulting  in  a  set  of  nmax  simultaneous  equations  in  the 

nmax  unknowns  Aq,  . . . ,  A^  max.    For  a  complete  shroud  we  must  evaluate 

qr     at  an  additional  point  to  give  us  nmax  equations, 
s 

The  solution  to  this  set  of  equations  provides  the  desired  representation  of 
the  circulation  distribution  and  permits  the  determination  of  the  side  thrust 
available  from  the  propeller  shroud  as  well  as  the  drag  induced  by  it. 

Weissinger  has  shown,in  Ref.  1,  and  it  is  shown  here  in  Appendix  C,  that 
this  method  yields  the  exact  value  of  total  lift  for  a  complete  ring  airfoil  in  the 
limit  of  very  high  or  very  low  aspect  ratios  (L/2ri/4)  and  yields  a  remarkably 
good  approximation  for  arbitrary  aspect  ratios  falling  between  these  limits. 

It  is  to  be  noted  that  this  method  gives  the  correct  solution  only  for  airfoils 
of  zero  thickness  as  it  takes  account  only  of  the  circulation  and  not  of  the  thick- 
ness.   The  effects  of  thickness  must  be  accounted  for  separately;  however,  their 
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effect  on  the  total  force  acting  on  the  surface  is  small  in  comparison  to  that  of 
circulation  as  shown  in  Fig.  22,  Ref.  2.    (Reproduced  here  as  Fig.  A-5.  ) 
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Sample  Calculations: 

Consider  a  propeller  with  the  following  characteristics: 

Dp  =  12.  0  '        SHP  =  6000 

(1-t)  =  0.80         (1-W)  =  0.70 

Vs  =  20  KTS  Va  =  Voo  =  14  KTS;  V^  =  23. 65  F.  P.  S. 

np  =  0.625         err  =  1.00 

N  =  200  RPM 

THP  =  SHP        nDerr  =  3750  HP  =  3~£- 
"    Xi  550 

CT  =  Thrust  Coeff.    =        0     — *- 

nR2lpV2 

Cm  =   2(THP)x550    .    p  =  1.99  siugs/ft3 
nRp  P(V/ 

c    _         2(3750)  x  550 7.  5  x  5.  5  x  106 

T  n(6.  0)2  (1.  99)  (23.  65)3  tt(3.  6)  (1.  99)  (2.  365)3  x  104 

Crp    =    1.  000 

Assume  the  circulation  distribution  over  the  propeller  is  constant;  radial  and 
axial  velocity  fields  for  this  case  are  tabulated  in  Tables  3  and  4,  Ref.  (8). 
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Shroud  Geometry. 

Consider  a  shroud  of  constant  chord,  camber  and  convergence  angle,  and 
assume,  for  the  purposes  of  this  calculation,  that  the  velocities  induced  by  the 
propeller,  at  the  shroud  vortex  sheet,  are  equivalent  to  those  induced  on  a 
cylinder  of  radius  Rp.    Take: 

%£1   =0.8;L  =  9.6/ 
Kp 

Choosing  NACA  a  =  0.  6  meanline  as  tabulated  in  Ref.  (7) 

c^  =  2.  58°   ;     c'|  x  =  y=  +  0. 1827 
£  |  x=   L  =  Q  0?37 

Referring  to  Table  3,  Ref.  (8)  at  £  ■  1. 0  ;  ^  ■  -0. 4  (shroud  leading  edge) 


^P        =  +  0.  08  ;  £      =  0. 08  (1. 385)  vo 


q      +v     =  1. 11  v 
^xp       °°  ° 


Referring  to  Table  4,  Ref.  (8): 

-♦ 

-2L    =  -0.086  ;    q_  =  -0.  086  (1.  385)  v 
v   CT  ^  °° 
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"*dVp 


•TRAILING-    VORTBTA 
SHEET 


■*-X 


FIGURE  3 


For  ideal  angle  of  attack: 


xp       °° 

O/p  =  -3  +  a4  =  6.12° 

-3  =  6. 12°  -  2.  58°  =  3.  54°  =  0.  620  Radians 

The  shroud  will  have  the  following  characteristics: 

rl/4  =  6'  °'  +  <°«  737>  <9-  6> +  °*  62  <9'  6/4>  =  6'  85?f 
L  =  9.60' 

®  =  90° 

Mean  line:   N.  A.  C.  A.      a  =  0.  6 
-p  =  3.  54° 

,         =W2_=    4^  =  0.70 
J/4       rj/4       6. 857 
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Boundary  Conditions  for  Determining  Thrust: 
(1)    at  9  =  ±  e;r(e)  =  0 


q       +  q  , 

(2)    at  X3/4;    J      -  ^    =  c      -  (-3)  =  0. 1827  -  0. 062 


From  Table  3,  Ref.  (8): 


-^L-  =  0. 155 


=  0.155  (1.385W    =  0.215v 

'  '      oo  oo 


From  Table  4,  Ref.  (8): 


-~E-  =  -0.136 

V       Crj, 
00       T 


q      =  -0. 136  (1.  385)  vo 


q  o  =  (0. 147  +  0.  188)  v     =  0.  335  v 
^rs  7co  c 


qrg  =  0.335(23.65)  =  7.  93F.P.S. 

Take: 

§!  =  0 

§2  =  tt/4 

l3  =  tt/3 
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The  set  of  simultaneous  equations  determining  the  four  coefficients  in  the 
cosine  series  representation  of  the  circulation  distribution  are: 

0.  025  A0  +  0.  0741  Ax  +  0.  0963  A2  +  0. 1739  A3  =  7.  93 
0.  0235  AQ  +  0.  0422  Ax  +  0.  0198  A^  -  0.  0450  A3  =  7.  93 
0.  0214  AQ  +  0.  0248  A1  +  0.  0165  ^  -  0. 153  A3  =  1.  93 
1. 0  A0        +  0  Al  -  1. 0  A2        +  0  A3  =0 

Solving  simultaneously. 

AQ  =  A2  =  338.  0 
Ax  =   -249.2 
A3  =  -84.6 

T/Q)  =  338  -  249.  2  cos  (0)  +  338  cos  (20)  -  84.  6  cos  (30) 

By  the  Kutta -Joukowski  Law 


'(aT) 


A  °°     %plx=oJ    „ 
-p £ r 


cos  0>p 


1/4 


,n/2 
-tt/2 


rfi  cos  0  d0 


ByRef.  (8): 


V 


V   CT 

OO         1 


=  0. 10  ;  q^  =  0. 1385  Vo 


x=o 


L(aT)   =  "188,500  Lbs. 


L,  Tv  is  positive  when  directed  radially  outward,  perpendicular  to  V       , 
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Determination  of  the  Induced  Drag: 

Di(aT)  =  -  T(0)  tan  o^ 

D^aT)  =  +  188,500  (0.045)  =  8,480  Lbs. 

The  induced  drag  is  in  a  direction  parallel  to  Vtotaj 


■*T(<» 


~P.*(o"> 


FIGURE  4 

-T(0)  =  -L/aTx  cos  Op  +  D.  (aT)  sin  aT 

-T(0)  =  188,500  (0.  9943)  +  8480  (0. 1067) 
-T(0)  =  187,  500  +  940  =  188, 440 
. -T^«  "L(aT)  =  188,  500  Lbs. 

-D^O)  =  "L/aTv  sin  (aT)  -  D.(aT)  cos  aT 

■Di(0)  =  20, 100  -  8,  400  =  11,700  Lbs. 


—  X 
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Form  Drag: 

For  an 

N.A. 

C. 

A. 

66 

-210  section 

ByRef.  (7): 

cD  © 

a  = 

0 

=  0, 

0094 

where 

CrJ 

Df 

jKV.  +  q^nr^L 


.*.  Df  =  0.  0094  (1.  0)  [<1.  215)  (23.  65)]2  tt  (6.  85  7)  (9.  6) 

Df  =  9.  4  (0.  826)  n  (6.  857)  (9.  6) 
Df  =  1600  Lbs. 

Taking  Df  as  acting  in  the  +x  direction 

Dtotal  =  Di  +  Df  ="11'  700  +16°0 
Dtotal  =  -9. 100  Lbs. 

Final  Predicted  Thrust  from  this  Shroud: 
Ahead  thrust  =  9,  000  Lbs. 
Side  thrust  =  188,000  Lbs. 
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Comparison  with  a  Conventional  Rudder: 
Assuming  that  a/6  -  5/7  where: 
a  =  angle  of  attack 
6  -  rudder  angle. 

For  a  rudder  of  Aspect  Ratio  =  3jfor  the  best  L/D  ratio  and  of  area  equal  to  the 
area  of  the  shroud: 

A  =  area  =  (6.  86  +  0.  48  )  tt  (9.  6)  =  221  Ft. 2 

AR  =  Aspect  Ratio  =  (2S)2/A  =  3 

S  =  Span 

S2  =  |  (221)  =  166  ;  S  =  12.  88' 


CL  = 

L 

cD 

D 

(p/2)  vf  A 

(P/2)  Vf 

A 

Figure  62,  Ref.  (9): 

6 

a 

6 -a 

Cl 

L 

QD_ 

D 

0 

0 

0 

0 

0 

0.01 

930 

22.4° 

16° 

6.4° 

0.830 

77,  000 

0.096 

8,900 

28.0° 

20° 

8.0° 

0.980 

91,000 

0.163 

15,120 
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FIGURE  5 


6  = 


TQ  =  Side  Thrust  =  L  cos  (6  -  a)  -  D  sin  (6  -  a) 
DQ  =  Drag  =  L  sin  (6  -  a)  +  D  cos  (6  -  a) 
22.4°: 


6  = 


6  =  0 


T0  = 

76, 000 

D0  = 

17,400 

8°:    ( 

Breakdo 

To  = 

88, 000 

i. 

27,  500 

1, 

To  = 

0 

Do  = 

930 
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SYMBOLS: 

An  =  Goefficient  in  series  representation  of  circulation  distribution 

C  =  Nondimensionalized  height  of  mean  line  above  the  nose  tail  line. 

C  =  Slope  of  the  mean  line  relative  to  the  nose  tail  line 

Cn  =  Drag  coefficient  =  — * — « *- 

"RP  2   PV 

Cl  =  Lift  coefficient  =  -    — ^— 


PV^  r1/4\ 


Cnm  =  Coefficient  in  series  representation  of  induced  radial  velocity. 

T 

Ct  =  Thrust  coefficient  = =-t ^ 

— ♦ 

D  =  Position  vector  from  vortex  element  to  pivitol  point. 

Df  =  Form  Drag 

Dj^a)  =  Induced  Drag 

E(cpp,k)  =  Incomplete  elliptic  integral  of  the  first  kind 

Fr  =  Radial  force  per  unit  arc  length 

Fx  =  Axial  force  per  unit  arc  length 

—»-»—♦ 

i,  j,k  =  Unit  vectors  in  Cartesian  coordinates 

Inm  =  Integral  occurring  in  the  series  representation  of  induced  radial 

velocity. 

9 

kA  =  Modulus  in  elliptic  functions 

K(k)  =  Complete  elliptic  integral  of  the  first  kind 

X    +  2 

t  =  — K =  constant  for  a  particular  pivotal  point 

L  =  Chord  length 

L(a)  =  Lift  force 

N  =  Propeller  R.  P.  M. 
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P  =  An  arbitrary  field  point  where  the  velocity  is  to  be  determined. 

q  =  Induced  velocity. 

— ♦ 

q^  =  Velocity  induced  by  the  bound  vortex  distribution. 

qr  ■  Radial  velocity  induced  by  the  propeller. 

qr  =  Radial  velocity  induced  by  the  shroud. 

qT  =  Velocity  induced  by  the  trailing  vortex  distribution. 

q  =  Axial  velocity  induced  by  the  propeller. 

q^s  =  Axial  velocity  induced  by  the  shroud. 

r  =  Radial  coordinate. 

r,cp,  k  =  Unit  vectors  in  polar  coordinates. 

R  =  Resultant  force  acting  on  the  shroud. 

SHP  =  Shaft  horsepower. 

t  =  Dummy  variable  =  cos  v  • 

t  =  Thrust  deduction. 

THP  =  Thrust  horsepower. 

T/g\  =  Side  thrust. 

v^  =  Speed  of  advance  through  medium. 

V  =  Free  stream  velocity. 

vl»v2»v3      =  Functions  used  in  evaluating  elliptic  integrals. 

— ♦  — »        — ♦ 

V„  =  V    +  q  _  =  axial  velocity  of  propeller  stream. 

op  °°     ^xp  J       r     r 

w  =  Wake  fraction. 

x,  §  =  Axial  coordinate. 

c^  =  Induced  angle  of  attack  at  the  1/4  chord  point. 

a0  =  Angle  of  attack  measured  from  incoming  stream  to  shroud  axis  in 

plane  of  symmetry. 

f3  =  Convergence  angle  of  mean  line. 

Y  =  Circulation  per  unit  length  (positive  in  the  positive  coordinate 

direction.) 
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F  =  Circulation  (positive  in  the  positive  coordinate  direction). 

e  =  — L-  =  Nondimensionalized  axial  coordinate. 
rl/4 

Err  =  Relative  rotative  efficiency. 

n^  =  Propeller  efficiency. 

0 ,  $  =  Polar  angle. 

0  =  Limiting  half  angle  of  shroud. 

X  -  x/rj/4  =  Nondimensionalized  axial  coordinate. 

v  =  Angular  variable  of  integration  =  rr  -  ($  -  6). 

§  =  Axial  coordinate. 

^pp,a,  ,k)  =  Incomplete  elliptic  integral  of  the  third  kind. 

p  =  Density  of  fluid. 

P  =  r/rw/  =  Nondimensionalized  radial  coordinate 

cpD  =  Amplitude  of  angular  argument  in  elliptic  functions. 

§  =  Angular  coordinate  of  the  pivotal  point. 


Subs 

cripts: 

B 

Due  to  bound  vortex 

n,m. 

■  P 

indices 

P 

Due  to  propeller 

r 

In  radial  direction 

s 

Due  to  shroud 

T 

Due  to  trailing  vortices 

X 

In  axial  direction 

e,cp 

In  tangential  direction 

1/4 

At  1  /4  chord  point 

3/4 

At  3/4  chord  point 
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APPENDIX  A 

DEVELOPMENT  OF  THE  INTEGRAL  EQUATION  FOR  THE 
VELOCITY  INDUCED  BY  THE  SHROUD 


I.    Shroud  Vortex  System: 


bound 

VORTEX 


CONTROL 
VOLUME. 


€ — C      cy 


TT*.A\UIN©    VOKTEX 


-X 


FIGURE  Al 


Consider  the  control  volume  enclosing  the  incremental  arc  of  the  bound 
vortex  distribution  lying  between  the  polar  angles  9  and  0  +  d0  of  Fig.  Al. 
Helmholtz*  vortex  theorem  states  that  vortex  filaments  cannot  end  in  the  fluid. 
This  requires  that  the  circulation  around  the  vortex  filaments  leaving  the 
control  volume  must  equal  that  of  the  vortex  filaments  entering.     This  theorem 
implies  conservation  of  vorticity  and,  in  terms  of  the  circulation  states  that: 


r(0  +  A9)  +  Yx(9)  r1/4  A0  =  T(0) 

.dT  (0) 
r(9)  +  C~di~0  ^  +  Yx(9)  rl/4  ^  =  r(9) 


Therefore: 
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A-l: 


W-^HP 


II.    Velocity  Induced  at  a  Field  Point  by  the  Shroud  Vortex  System: 

By  the  law  of  "Biot-Savart": 
A-2: 

dq  =     1      r(ds  XD) 

W         |D|3 

where:       dq  =  incremental  velocity  induced  at  a  particular  field  point  or 
"pivotal  point", 
as  =  directed  incremental  element  of  the  vortex  line. 

T  =  magnitude  of  the  circulation  around  the  vortex  element  3s. 

— » 

D  =  vector  displacement  from  the  vortex  element  to  the 

pivotal  point. 

IntegratingA2,   it  follows  that: 
A-3: 

%<*••*-£  J8       r(9)ri/4^3X5 
^  4tt  J-@  |D|3 

and 
A-4: 

®     co 


qT(r,,,X)  =  ^  J  J       Yx(9)r1/4(3xXD)de 


■®  o  iD|3 


If  the  pivotal  point  coincides  with  a  point  on  the  vortex  system,  one  or  both  of 
the  integrals  must  be  evaluated  as  the  Cauchy  principle  value  to  exclude  the 
singularity  at  the  pivotal  point. 
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III.    Integral  Equations  for  Velocity  Induced  by  the  Bound  Vortex  Distribution: 


■*-X 


FIGURE  A2 


"P"  is  an  arbitrary  pivotal  point  where  the  induced  velocity  is  to  be  deter- 
mined,   p  has  cartesian  coordinates  (x,  y,  z)  and  polar  coordinates  (r,  $,x). 
By  the  law  of  Biot-Savart: 

A -4: 

deXD 


3%  =  4^  r  rl/4 


D 
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In  cartesian  coordinates: 
A-5: 

d9  =  d9  [Of-  sin  9j*+  cos  0k] 

r  =  r  [Oi  4-  cos  $  j  +  sin  $  k] 
A-6; 

D  =  (x)~T  +  [r  cos  $  -  ri/4  cos  0]j  +  [r  sin  $  -  ri/4  sin  8]k 

30XD  =  -d0    sin  0  (r  sin  $  -  r^M  sin  9)  +  cos  9  (r  cos  $  -  rj/4  cos  0)ii 
+d9  [x  cos  0]"j 
+d9  [x  sin  0]  k 

|D|  =  [x2  +  r2  +  r2l/4  -  2r  r1/4  cos  ($  -0)]1/2 

Nondimensionalizing: 
A-7:      ' 

r  -      o     •  X  -     1 

-  P  >  zr~r  -  X 


rl/4  rl/4 

A-8: 

|D|   =  r1/4[X2+p2  +  l  -2pcos($-0)]1/2 

A-9: 

cfeXD  =   r!/4  |-[p  cos  ($  -0)  -  l]d9  i+  X  cos0d0"J  +  X  sin0  d0ki 

Substituting  in  A-7,  A-8,  A-9  into  A-4: 


da       =  -      H9)  [pcos($-9)-l]de        _        7 

%X"      4nrl/4       [X2  +  p2  +  l-2pcos(*-0)f/2 


dj  r(8 )  X  cos  0  d0 


by      4nr1/4  [X2  +  p2  +  1  -  2p  cos  (§  -  0)]3/2 
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dq 


B  T(8)X  sine  d9 

4nr1/4  [X2+p2  +  l  -2pcos  ($-9)]3/2 


dqbr=dqb  "tt 


dqbrr  = 


T(9 )  X  cos  9  cos  $  d9 

4nr1/4  [X2  +  p2  +  l  -2pcos  ($-9)]3/2 


T(9)  X  sin  9  sin  $  d9 


+  —    w/      ^  u  ^      w 

4nrl/4  C^  +  P24"1  "2pcos  ($-9)] 


372" 


dq 


Xr 


br      4nr 


1/4 


r(9)cos  ($-9)d8 
[X2  +  p2  +  1  -  2p  cos  (f  -9)]3/2 


Integrating  over  the  shroud: 


A-10: 


— ♦ 

X? 

qbr  " 

4nr1/4 

Similarly: 

A-ll: 

I 


© 


r(9)  cos($  - 9)d9 


!/4     J_@     [x2+p2  +  l-2pcos  ($-9)] 


3/2 


0 


-l 


r(9)[pcos  ($-9)  -  l]d9 
*bx  -  4nr1/4    J_@   [x2+p2  +  1.2pcos(^e)]3/2 


-r(9 )  X  cos  9  sin  $  d9 


M         [4nr1/4  [X2+  p2  +  l  -2pcos  ($-9)]3/2 


+  T(9)  X  sin  9  cos  cp  d9 


4nr1/4[X2  +  p2  +  1  -  2p  cos  ($  -  6)]  ' 


*  "P 
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A-12: 


*b§       4nr 


1/4 


8 


-X-.cp       P      r(8)sin_($-_9),d9. 


-©  [X2  +  p2  + 1  -  2  p  cos  (§  -  0 )]  ^ 


Note:    q       is  of  little  interest  in  this  problem  and  is  included  only  for  the  sake 
of  completeness. 

IV.    Integral  Equations  for  the  Velocity  Distribution  Induced  by  the  Trailing 
Vortex  Distribution: 


**- 


« — fr— £ — £ — Q 


c3p  \^ 


FIGURE  A3 
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ByA-1: 

By  the  law  of  Biot-Savart: 


a-       \  f    '  \     dr(e)  ^         .dfxB 


A-13: 


dq*  =  -  x  rageTu odgx5 

dqT  4n     V  de     ya9      iD|3 


D  =  (x  -  §)  i  +  (r  cos  $  ~  *  1  /4  cos  9 )  j 
+  (  r  sin  $  -  rw,  sin  e)k 

— »  — ♦ 

d?  =  d§  i 

d?  X  D*  =  (0)f  -  d§  (r  sin  $  -  rA  /4  sin  9)7 
+  d§  (r  cos  $  -  rw4  cos  9)k 

| D|   =  [(x  -  ?)2  +  r2  +  r2/4  -  2r  r1/4  cos  ($  -  9)] 

Nondimens  iona  liz  ing: 
By  A- 7: 

r         „  .     x         , 

=  p  ;  - —  =  a. 


rl/4  rl/4 

A-14: 


Define      § 

— —  =  e 

rl/4 
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A-15; 

|D|  =  r1/4[(\-e)2  +  p2  +  l  -2pcos($-e)]1/2 
A-16: 

d§XD   =  -  rwJ  de(p  sin  $  -  sin  0)"J  -  de  (p  cos  $  -  cos  0)k  J 

Substituting  A-14,  A-15,  and  A-16  into  A-13: 
dqrpx  =  0 

2n   v  =        I        ^dr(8)^        (psin$  -  sine)  ded9? 

^  "  4nr1/4  k   dG    ;[(,.e)2+p2  +  1.2pcos($.e)]3/2 


-i  1         ^dr(e)  \  (p  cos  $  -  cos  e)dede  k 

<%*  -    -4nr1/4    V.   d8    J  [(,  .e)2  +  p2  +  \  .  2p  CQS  (f  .e)]3/2 


dqTr  =  dqT  •  -~  =  dqT  •   [cos  §  j*  +  sin  $  k] 


ri       -♦  _  1_      dr(e)  [+sin  ($-8)]ded9  r 

^Trr-4nr1/4       de      [(,  _  e)2  +  p2  +  x  .  2p  cos  ($  _0)]3/2 


dqT$  =  dqT  •7^7  =  dqT  •   [-sin  $7  +  cos  $  k] 
l<*Pl 

,       -  1_     dT(9)        [P  -  cos  ($-9)]  de  d9  cp 

aW-- 41^/4     d9       [(,.e)+p2  +  1.2pcos(,.e)]3/2 

integrating  over  the  shroud: 

A- 17: 

— ♦  @  00 

*TY  "  4nrr  ,  (6)  sin($-9)ded9 

^       4nil/4    "U    ""a^[fc-e?  +  ^  +  l-2pcos<*-e>]4/2 
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A-18: 


TS 


@      00 

■  .-? r  r 

4TTri/4-iJ0 


dr 


(9)     .  .  .  [p  ~  cos  ($  -  0)]  de  de 

de 


B\-e)2+  p2  +  l-2pcos  (S-e)]3^ 


A-19: 


q-px 


=  0 


Considering  A-17: 


lTr      4nr 


©  dr 


1/4    J-(B) 


de 


fil.i 


sin 


C*-e){ 


de 
[(\-e)2+k2]3/2 


}de 


where  k     =  £r  +  1  -  2p  cos  ($  -6)  can  be  treated  as  constant  when  integrating 
the  integral  in  brackets  (k    >.  0).    Integrating  the  bracketed  integral  by  parts: 
Let  u  =  \  -  e;du  =  -de 
LIMITS:    e  =  0  =>  \x  =  \ 

£    =    »=)     U=-oo 


J 


de 


[(\-e)2  +  k2)3/2 


du 


[u2  +  k2]3/2   ~kV+k2)1/2 


kV+kV/2  +k2 


A-20: 


rA 


„®  dr 


irTr      4nr 


(9) 


1/4    J-© 


de 


sin  ($  -  e)  de 

[p2  +  l  -2pcos  (§-e)]  [X2+  P2  +  1  -2pcos($-6)]1/2 


r      r@dr(e) 

4nrl/4  J-@     d0       [p2  +1  -2pcos($-e)] 


sin  (§  -  9)d0 
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V.     Complete  Integral  Equations  for  the  Velocity  Induced  by  the  Shroud: 
Combining  A- 10  and  A-20: 


^s  =  qbr  +  qTr 


A-21: 


\r 


© 


r(e)  cos  ($-e)  de 


rs       4nr1/4     J_@    [x2  +  p2  +  x  -  2p  cos  ($ -0)]3/2 


X@    dr0) 


4nr- 


1/4    J. 


sin  (cp  -  8)  d9 


© 


d9 


[p2*!  -2pcos($-9)]  [X2+p2  +  l  -2pcos($-9)]1/2 


+■ 


£®  dr(9)  sin($-9)d9 


4nr 


r       I  (9)  sin  (3 

1/4    *-©    de       [p2  +  l- 


2p  cos  ($-0)] 


ByA-11: 


Scs  "  qbx  " 


8 


r(e)[pcos  ($-0)  -  i]d0 

4TTrl/4     J-@    [\2  +  p2  +  1  -2pcos($-0)f/2 
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APPENDDC  B 

SOLUTION  OF  THE  INTEGRAL  EQUATIONS  FOR  A  PIVOTAL 
POINT  COINCIDENT  WITH  THE  TRAILING  VORTEX  SHEET 

I.    Reduction  of  Eg.   (A  -21)  to    Inegrable  Form 

For  a  pivotal  point  on  the  trailing  vortex  sheet,  r  =  r^  m,   p  =  1,  and  the 
integrals  involving  dF/Q\/d0  which  arise  from  the  trailing  vortices  must  be 
evaluated  as  Cauchy  Principle  Values. 
A -22: 

,®        r,0*  cos  (§-e)de 


<lr 


(i) 


3/2 


s       4nr1/4     J_@    [y2  +  2-2  cos  ($  -  9)] 


® 

t    Xr       ft    dr(ej 


8TTrl/4    ^     de        [l-cos($-e)][?,2  +  2-2cos($-e)]1/2 


-»  ®  dr 

+       r         L      ai(9)        sin($-9)d9 

8nrl/4   i®     de       [l-cos(§-9)] 


Assuming  that  the  circulation  distribution  can  be  represented  as  a  com- 
plete Fourier  series: 


r(9)  =   Y    An  cos  (n9)  +  Bnsin  (n9) 


n=o 

The  boundary  conditions  on  this  representation  are: 

a.  F(9)  is  an  even  function  of  9  for  a  shroud  symmetrical  about  the 
9=0  plane. 

b.  r(@)  =  r(-®)  =  0  (®  ?  n);  that  is  the  circulation  goes  to  zero  at  the 
tips  of  the  shroud. 

By  (a)  Bn  =  0  for  all  n. 
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A-23: 

03 

T(e)  =    £  A^cos  (n0) 


n=o 


=  -    V  nAnsin(n9) 


de 

n=o 


A-24: 

Define   -     '  -  =  I 


\2  +  2 


Substituting  A-23  and  A-24  into  A-22  and  summing  over  the  integrals: 

A-25: 

r  X  1  An        e 

n=o  f      cos  (n0)  cos  (g  -8)d9 


a       _  n=o  f      cos  (n8)  cos  ($  -9)d( 

rs      8/2  n  r1/4      {@    It  -  cos  (*  -  9)]3/2 


?XI^         © 


n=o  £         sin  (n9)  sin  ($  -  9)  d9 

&/2nr,  ,„      .p)   [l-cos(§-e)][<t-cos(§-0)] 


1/4     -® 


172 


?InAn 


£    sin(n9)gfii^ 
«T  [1  -cos($  -£ 


8nr1/4       le  Cl-co8(#-9)] 

Transformation  of  variables: 

n9  =  n$  -  n($  -  9) 

cos  n9  =  cos  n($  -9)cos  ncp  +  sin  n($  -9)  sin  ncp 

sin  n9  =  cos  n($  -  9 )  sin  ncp  -  sin  n  ($  -  9 )  cos  ncp 
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Define: 

vnn-^-e);  dv=d9 
sin  (<£  -0)  =  sin  v 
cos  ($  -0)  =  -  cos  v 

cos  n0  =  cos  n(n-  v)  cos  n$  +  sin  n(n-  v)sin  n$ 
sin  n0  =  cos  n(n  -  v)  sin  n§  -  sin  n(rr  -  v)  cos  n$ 
since 

n(n  -  v)  =  nn  -  nv  ; 

cos  (nn)  =  (-l?1 

sin  (nn)  =  (0)™ 

cos  n(n-  v)  =  (-l)n  cos  nv  +  (0)n  sin  nv  =  (~l)n  cos  nv 

sin  n(rr  -  v)  =  (Op  cos  nv  -  (-l)n  sin  nv  =  -(-l)n  sin  nv 

.'.cos  n0  =  (-If1  [cos  nv  cos  n§  -  sin  nv  sin  n$] 

sin  n0  =  (-If1  [cos  nv  sin  n$  +  sin  nv  cos  n§] 

Revised  limits  of  integration: 
when 

0  =  +©  ;  v  =  tt  -  (§  -  ®)  =  n+®  -  $  =  v2 
when 

0    =    -©;    V=TT-($+©)  =  TT-®-§    =  Vi 

Note:    Restriction:    vj»  \£  7*  tt     §^±© 
Further  defining: 

c       =  /  1>n+l  X  C0S  (n^)An  -      u   f  ^  cos  (nv)  cos  v  dv 

nl  8/2  tt  r1/4  ™  "  Jvx         [£+cosv]3/2 


•38- 


Cn2  -  ("I)" 


X  sin  (n$)A 


8/2nr 


1/4 


T       _   r^  sin  (nv)  cos  v  dv 

^  =  J       771         iVi 

vi     11  +  cos  v] 


Cn3  s  ("« 


.,  \  n  sin  (n$)  A                        „yo          ,     »     .        , 
n+1  x     '     n  T  v   z  cos  (nv)  sin  v  dy 

n3  "  J       7T.  777         7172 


I 


8/2iTr1/4 


v^    [1+cos  v][t+cos] 


cn4  -  (-D 


n+1  X  n  cos  (n$)  An 


8/2  nr 


1/4 


T  f^        sin  (nv)  sin  v  dy 

\>l    [1  +  cos  yj[t+cos  v]  ' 


Cn5  =  (-1) 


n+l  n  sin  <n*>An 


8nr 


1/4 


"nS" 


vl 


^2    cos  (ny)  sin  y  dy 
[1+cos  v] 


C       =  (-1)114"1  n  C0S  ^n$^An 
8TTrl/4 


^6  " 


^2  sin  (ny)  sin  v  dy 
[1+cos  y] 


vl 


Based  on  the  above  definitions  A -2  5  reduces  to: 


A-26: 


^s  =  ?  I     I      Cnm  W 
n=o  m=l 

Solutions  to  the  above  integrals  are  given  in  the  following  equations:  (nj£  3) 


-39- 


Integral 

E^ 

Integral 

Eq. 

m=  1: 

% 

A-43 

m 

-4: 

^4 

A-47 

111 

A-44 

114 

A-48 

hi 

A-45 

x24 

A-49 

hi 

A-46 

134 

A-50 

m  =  2: 

!02 

A-27 

m 

=  5: 

ks 

A-35 

Il2 

A-28 

hs 

A-36 

*22 

A-29 

hs 

A-37 

*32 

A-30 

135 

A-38 

m  =  3: 

*03 

A-31 

m 

=  6: 

X06 

A-39 

*13 

A-32 

*16 

A-40 

X23 

A-33 

he 

A-41 

Z33 

A-34 

he 

A -42 

-40- 


Note:     In  evaluating  the  integrals  lni  and  1^  which  involve  elliptic  functions 
the  following  relations  are  used: 


k2  =      2       ,  k/2  _  ,.2  =  ill 


2  g        /ttt       ,    ,  ,  '2       21 


e-  =      z      •  _£.  =  /I+T  •   i  +  k     =   l^ 


/(4+l)(l+COS  Vp) 

cpp  =  am  up  =  sin       ^      2(<t  +  COSVp) 

Note:    See  Figure  (  A -5)  for  quadrant  selection. 
sn(up,  k)  =  sin  cpp  =  sn  up  ;  en  up  =  cos  cpp 
dn  u  =  J  1  -  k2  sin   cp 


1  «.  sn  u  ,  sn  u 

ns  u  =  tn  u  =  sc  u  =  sd  u  =  +3 

sn  u  en  u  dn  u 


1  en  u  ,  en  u 

nc  u  =  cs  u  =  — ■ —  cd  u  = 


tn  u  =  sc  u  = 

sn  u 
en  u 

en  u 

cs  u  =  

sn  u 

,            dn  u 
ds  u  =  

en  u  '  sn  u  dn  u 


,1                         ,            dn  u  ,  dnu 

nd  u  =  -3 ds  u  =  dc  u  = 


dn  u  sn  u  en  u 


VT(u)  =  -i-  [E(u)  -  k2  sn  u  cd  u] 
k'z 


V,(u)  =  — K-   [2(2-k2)E(u)  -k/2u  -k2snucdu(k/2nd2u  +  4-2k2)] 
3k'* 


V  (u)  =  4(2  -k2)  V2  -  3Vi  -  k2  snu  cdu  nd4u 
3X1    '  5k/2 


-41 


Summary  of  Integral  Values: 


IQ1  =_!_  |U1  -  u2  --Lj.  [E(ux)  -  E(u2)]} 

A+l  * 


(A-43) 


In  =  ^kill  ^(cp^k2^)  -  n(cp2,k2,k)l  -- 2i-  (Ul  -  u2)  -  tClQi]  (A -44) 

A+l     L  J    /t+i 


I21  =  2t2(Ul  -  u2)  +  4^(1 -t)[Vi(u1)-V1(u2)]  + 2(1- 1)2  [V2(Ul)-v2(u2)] 

-2^[Iu]-[loi]  (A -45) 

I31  =  4t3(u1-u2)+12t2(l-^)[V1(u1)- V1(u2)]+ 12^(1 -t)2[V2(u1)-V2(u2)] 

+4(l-^)3[V3(u1)-V3(u2)] 
-2^[I2l]  -3  [In]  -2£[loi]  (A-46) 


102  =  0 


I12  =  2[a+tP  +  Ul+t)  2  ] 


cos  vi 


(A-27) 
(A-28) 


cos  V2 


^3/2 


l22  =  4{^i^+4[Mtpi  .^+t)ii} 

La+t)5      L      d  JJ 


cos 


COS 


vl 


^ 


(A-29) 


I32  »  {-^  (3t2  -  4tt  +  SI2  -  |^)  Ol+t)2  -  8(t3  +  2*,)  (t+t)"*} 


u    cos  vi 


cos  V2 


(A-30) 


L03 


A-i  l    L 


/Z+t  -/FT  \  lcos  V1 

/t+t  +  /rnr 


} 


cos 


^2 


(A-31) 


-42- 


cos 


vl 


I13  =  2  /{,+t 


-[l03l 


cos  y2 


(A-32) 


%  =  (4-^-4)^ 


cos 


COS 


vi 


^2 


+  Cl03l 


(A-33) 


I33  =  "jg-    3t2  -  5t  -  4<tt  +  8t2  +  lOt  +  —■     /I+T 


h4  =  0 


cos 


"cos 


vl 


■CI03] 


^2 


I34  =  8^"  I)2  /t+r{v!(u2)  -  V^u^  +  ^j-  [v2(ui)  -  V2(u2)] 

+  ZTT  [V3(u2>  "  v3<ul)]}  "  2<^+1)  CI24]  "  (4t+l)  [I14] 


I05  =  tn 


pl+COS    Vn      -1 

_l  +  COS   \^    - 


I15  =  cos  vi  -  cos  \>2  -  [I05] 


I25  =  cos    V!  -  cos    \>2  -  2[Ii5]  "  [I05] 


I35  =  §  (cos3  V1  -  cos3  V2)  -  2[I25]  -8[I15]  -4[Iq5] 


(A-34) 


(A-47) 


r  9  ~i 

I14  =  2/1+I  I  ui  -u2  +  E(u2)  -  E(ui)+y£j  (snui  cduj  -snu2  cdu2)|         (A -48) 

L>4  =  I  (f  1)  /£+T  {2u2  -  2ux  +  ^  [E(ui)  -  E(u2)] 

+irr  [Cnd2ui  "!^r)snui  cdui  -(nd2u2  -$0* 

snu2  cdu2    ]•  -  2  [I14] 


(A -49) 


(A-50) 

(A-35) 
(A-36) 
(A-37) 
(A-38) 


-43- 

I06  =  °  (A-39) 

Il5  -  sin  vi  "  sin  V2  -  vi  +  V2  (A-40) 

*26  =  ^  sin  ^2  ~  ^  s*n  vl  +  2  (sin  2vi  "  sin  2y2)  -  (±  v^)  +  (±Vo)  (A-41) 

Note:      Sign  is  negative  if  vp  lies  in  the  first  or  second  quadrant  and  positive 
for  third  and  fourth  quadrants. 

I36  =  |  (sin3  V2  -  sin3  Vj)  -  2[I26]  -  [I16]  (A-42) 


-44- 


U.     SOLUTION  OF  THE  GEOMETRIC  INTEGRALS,  lnm  IN  EQUATION  A-26 

Integrals  L^  and  L^  are  elliptic  integrals  and  are  most  readily  solved  in 
terms  of  Jacobian  Elliptic  Functions.    We  will  defer  this  step  until  last  and  first 
consider  integrals  In2»  In3»  In5»  and  In^  which  can  be  reduced  to  elementary 
forms. 

Using  the  algebraic  transformation: 

cos  v  =  t  ;  sin  v  =  /l  - 12 
-sin  v  dy  =  dt 


dt 


A- 

■  t2 

n 

cos  ny  =  gn(t) 

0 

1 

1 

t 

2 

2t2-l 

3 

4t3  -  3t 

sin  ny  =  fn(t) 

0 

/(1-t 

2) 

2t/(l 

-A 

3/(1 

-t2) 

-4(1- 

t2)Jv~- 

t2) 

=  (4t2 

-dAT- 

•t2) 

Considering  L^: 


Ln2 


v2 

sin  n  vcos  yd  y 

Jvl  [Ucosv]3/2 


1 

=  1 


COS  V 


y  (t)  t  dt 

n 


cos  \^    (l-t2)1/2  U+t]3/2 


(Note  reversal 
of  limit  with 
sign  change) 


•45- 


For  n  =  0:    fR(t)  =  0 
A^7:    I02  =  0 


For  n  =  1:    ijt)  =  /Tt2" 


[12=+J 


cos  vi  tdt 


cos  V2      [£+0 


372 


By  119,  Ref.  (3): 


x  (a+bx)  dx  =  -7T 

b2    L 


(aW4±n)/2       a(a+bx)<2±n)/2 


4±n 


2±n 


with  a  =  £,  b  =  1,  ±n  =  -3. 


A-28:    I12  =  2[(^  +  t)1/2  +  -C(t  +  t)"1/2j 


COS  V 


cos  V2 


For  N  =  2:    fn(t)  =  2t  /(1-t2) 


X22  "2 


-  cos  vi      t2  dt 


cos  v2    [£+t] 


3/2 


Integrating  by  parts  with  u  =  t  ,  du  =  2t  dt 


dv  = 


dt 


(l+t) 


3/2 


v  = 


-2 


a+o 


T/2 


*22  "  2 


2f 


a+o 


172 


cos  V} 


cos  Vi 


+  4 


t  dt 


cos  V2 


cos  V2    (-t+t) 


172 


By  119,  Ref.  (3):    with  a  =  I,  b  =  1,  ±n  =  -1 


-46- 


A-29:     I22  =  4  < 


-t 


(t+t) 


172 


+  4 


U+t) 


3/2 


1(1+ 1) 


1/2 


cos  vi 
cos  V2 


For  n  =  3:    fn(t)  =  (4t2  -  1)  /l-t2 


%_-* 


COS  VI  o 

tddt 


372 


cos  vi 


t  dt 


cos  v2    [£+_]'  cos  v2  [£+t] 


J72 


o  9 

Integrating  the  first  by  parts  with:     u  =  t  ,  du  =  3t    dt 


dv  = 


dt 


(t+t) 


372" 


a+t) 


172 


The  first  integral  equals 


I 


2V 


(l  +  t) 


172 


cos  vi  cos  vj      t2  ± 

+  o 


COS  V2 


cos  V2 


(-L+t)1/2 


> 


By  112,  Ref.   (3): 


r   x2 dx         2(8a2  -  4  abx  +  3b2x2     ,-— - 
_____  =  __ /a  + 

/a+bx  15  b0 


/  a  +  bx 


with  a  =  t,  b  =  1: 


8tv 


j-72  +  if   (Sl2  '  4U  +  3t2)  (^+t)1/2 


<*+t) 


The  second  integral  is  just  A-28  .*.  combining: 


COS  V} 

cos  V2 


A  -30: 


L32 


=  <^  (312  -4^  +  8-t2  -|)(t+t)1/2  -  (8t3  +  2t)a+t)'1/2 


COS  Vj 

COSV2 


•47- 


Considering  1^: 


Ln3 


*        cos  nvsinv  dv 


yi      [1  +  cos  v]  [l+  COS  vj 


72 


COS  v 


gn(0  dt 


cos  v2     [1+t]  [<t  +  ] 


172 


For  n  =  0:    gn(t)  =  1 


L03 


cos  v. 


dt 


Jcos  v2      [l+t][£  +  t] 


172 


By  56,  Ref.  (3): 


I 


dx 


1  b'/u  -  A5' 


v/u        /lEb'        b'/u  + /IE' 


with 


v  =  1  +  t  =  a'  +  b't  ;     a'=  lb'=  1 
u  =  £  +  t=a  +  bt;  a  =  £,  b  =  1 

k  =  ab'  -  a  'b  =  (l  -  1) 


A-31: 


Ir>3  - 


03  -  /rrr    i tn 


/CTt  +  /-t  -  l 


cos  vi 
cos    2 


For  n  =  1:    gn(t)  =  t 


COS    Vn 

1  =  r  tdt 

4osv2    [l+t][<L  +  t]1/2 


■48- 


Using  v  and  u  as  defined  in  IQ3,  t  =  v  -  1 


L13 


- 


cos  Vl  cos  vl  cos  VI 

(v-l)dt         r  dt 


cos  V2      v/~u" 


j 


cos 


dt 


V2     /u~       cos  V2     v  /*u" 


=  52/1 


+  t  - 


/T"! 


£n 


/I+t     "    Orl 


/Z+t  +  /t-i 


cos  vi 
cos  V2 


A -32: 


I13  =  2/I+t 


cos  v^ 
cos  V2 


L03 


For  n  =  2:     gn(t)  =  2t2  -  1  =  2(v  -  l)2 
gn(t)  =  2v2  -  4v  +  1 


-COS  VI    „    ,. 

cos  v2    /  u 


rcos  Vl    ^        rcos  Vl      dt 


COS  V 


:os  vo    v/^7 


since  u=<l  +  t=l  +  t-(l-<t)  =  v-(l-<6) 
dt  =  dv 

1  +  COS  Vi  j 

J-        v  dv 


the  first  integral  =  2 
By  111,  Ref.  (3): 


1  +  cos  v2    A^-  l)+v 


2  r_i^L-=  -2<2a:bx)/TTb7 

J  /a+bx  3b 


with:  x  =  v,    a  =  (£  -  1),    b=l 

the  first  integral  =  -  £  (2£  ■-.  2  -  v)  ^(t-l)  +  v 


1+cos  vl 
1  +  cos  V2 


-49- 


reverting  back  to  t: 

4 


=>  -  |  (21  -  3  - 1)  /t+T 


COS   V} 

cos  \>2 


the  second  integral  =  -8/1  +  t 


COS  Vi 

cos  \»2 


the  third  integral,  from  A 


30 


T-i 


In 


Jl+t  -H-\ 


Ji+t +/1-1 


cos  \>i 
JJ  cos  \>2 


Combining: 


hs-  (5-¥-*y*Z*£=i 


I 


In 


Jl+t 

- 1/4,-1 

yi+t 

+  A-1 

J 

cos  vi 
cos  V2 


A-33: 


cos  vi       r       -. 
[03J 


cos  V2 

For  n  =  3:     gn(t)  =  4t3  -  3t 

=  4(v-l)3  -  3(v-l) 
gn(t)  =  4[v3  -  3v2  +  3v  -  1]  -  3v  +  3 
=  4v3  -  12v2  +  9v  -  1 

By  the  procedure  followed  in  solving  I23: 

1  +  cosv!       y2dv 


"as-"4 


1  +  C0SV1  vdv 

*!l  +  cos  V2    */(<&- l)  +  v  Jl+cos  V2    </(£-!)  + v 


=  - 12  J 

+  v  Jl+( 


+  9 


cos  V!    dt         pcos  vx    dt 


cos  v2    Ai        cos  V2   v*/u 


-50- 


By  112,  Ref.  (3): 


*    x2dx      .  2(8a2  -  4  abx  +  3b2x2) 


•  a  +  bx  15b 

with  x  =  v,  a  =  £  -  1,  b  =  l,  the  first  integral  yields 


/a  4-  bx 


8[8«--l)2-4(Y5Pv+3v2   ]/(I7Ty?7 


1  +  cos  vi 
1  +  cos  \»2 


reverting  back  to  t: 


=>  jg    SI2  -  161  +  8  -  41,  +  4  -  4lt  +  4t  +  3  +  3t2  +  6t]  /ZTt 


COS  Vi 

cos  V2 


_8_ 
15 


3f 


+  lOt  -  4  It  +  8-L2  -  201  +  15  J  /Z+t 


COS  VJ 

cos  V2 


The  second  integral  yields,  from  the  development  of  A-33: 
8(21,  -  3  -  t)  ST+i 


cos  vi 
cos  ^2 


The  third  integral  yields: 


18/T+t 


COS    V} 

cos  V2 


The  fourth  integral,  from  A-30  yields: 


1       .     A+t  -Vl,- 1 
Jl-l         Ji+t  +  ji-i 

cos  vi 
cos  V2 

combining: 

■51- 


L33 


il  [* 


+  lOt  -  15t  -  4  £t  +  8V-  -  20*1  +  30t  +  15  -  45  + 


¥\ 


Jl  +  t  - 


J I  -  1 


-In 


Jl+t   -  /t-1 

/Z+t  + /t- 1 


cos  v 


cos  v. 


=  < 


15 


3t2  -  5t  -  4  £t  +  8t2  +  104,  +  ■=£ 

4 


/Z+t 


A-i 


-tn 


/I+t  -  /£-i 
_/Z+t  +  /£- 1_ 


i 


cos  \>i 
cos  V2 


A-34: 


I33  =  ^|  (j3t2  -  5t  -  4  -tt  +  8 12  +  10 £  +  ^  ^)/Z+7 


cos  v 


1 


cos  V2 


-  [103] 


Considering  In5: 


I 


n5 


^2 


cos 


Vi 


cos  nv  sin  y  dy 
[1+cos  y] 


cos 


Vl     gn(t)dt 


^2 


[1+t] 


For  n  =  0:    gn(t)  =  1 

cos  y^        ^t 


105  = 


cos  y2 


[1+t] 


=  to(l+t) 


COS  yi 
cos  \>2 


A -35: 


I05  =  to 


1+cos  v\ 
1  +  cos  \>2 


For  n  =  1:    gn(t)  =  t 


-52- 


L15 


rcos  vi   t  dt 


cos 


v2 


[1+t] 


By  32     Ref.   (3) 


I 


x  dx  1 

(a+bx)    =  .2  L 


a+bx  -a.  In. 


(a+bx)] 


with  x  =  t,    a  =  b  =  l 


i15  =  [i  +  t-ta(i+t)]|cosvi 

I  cos  V2 


*15  =  cos  vl  "  cos  ^2  "  ^n 


1  +  cos  vi 


1  +  cos 


V2 


A -36: 


I15  =  cos  vi  -  cos  V2  -  [l05l 


-  o<-2 


For  n  =  2:    gn(t)  =  2tz  -  1 


*25  =  2 


COS  Vi  cos  vi 

t2  dt  r  dt 


cos 


^2 


:z  dt  [* 

:i+t]  "J 


cos  y2 


[1+t] 


By  36,  Ref.   (3): 


J  f+ff  =  4  [ I  (a+bx>2  "  2a  <a+bx>  +  a2  <n  (a+bx)] 


with  x  =  t,    a  =  b  =  1,  the  first  integral  of  I25  =) 


2[(1+2t)      -2(l+t)+<n(l+t)] 


cos 


vl 


cos  V2 


-53- 


combining  with  -Iq5 


I25  =  2{4--t-|  +  ita(l+t)} 


cos  vi 
cos  V2 


1+COS   Vi 


I25  =  (cos2  Vj  -  cos2  v2)  -  2  (cos  V!  -  cos  v2)+  *n  [1+cos    *  ] 


V2 


A-37: 


I25  =  cos2  Vi  -  cos    \>2  -  2CI15]  -  [I05] 
For  n  =  3:    gn(t)  =  4t3  -  3t 


L35 


-J 


cosvi    t3dt 


cos 


VI 


cos  v2    t-l+t] 


t  dt 


cos 


v2 


Cl+t] 


By  the  recurrence  relation,  Equation  2-1112,  Ref.  (4): 


1 


xn  dx 


.-n 


a   x11"1  dx 


(a  +  bx)     "  nh       b  (a  +  bx) 


with  x  =  t,  n=3,  a  =  b  =  1  the  first  integral  of  I53  yields: 


t2dx 
3    "  J  IT+tT 


From  I25  and  1^: 


I35  =  4(4-[li^-   -2(l+t)  +  to(l+t)]} 


COS  Vi 


COS 


V2 


-  3 


1  +  t  -  £n(l  +  t) 


cos 


cos 


^2 


V2 


-54- 


combining: 


I35  =  3  [(cos  vi )3  -  (cos  V2)3  J  -2 1  cos2  vj[  -  cos2v2j 


1  +  cos 


i-l  -t-  cos  vj  -1 
+  (cos  Vl)  -  (cos  v2Hn[1  +  cQs        I 


A -38: 


I35  =  I  [(cos  V2)3  -2[l25]  -3[l15]  -4[lo5[ 


Considering  1^: 


v2 
T  f       sin  (ny)  sin  v  dy 

n6  =  J  LI  +cos  y] 

vi 

For  n  =  0:    fn(t)  =  0 

A-39: 


cos 


Vl  f  n(t)  dt 


cos 


V2 


[1+t] 


106  ■  0 
For  n  =  1:    fn(t)  =  A-tz 


I 


16 


cos 


COS 


vl 

^2 


COS   Vl     f  , 


rrtTdt  =  J Vi 


[i+t] 


cos  v2 


+t 


dt 


By  60,  Ref.   (3): 
'/v 


r/v  1    _       k 


J^rtan-1 


T7^ 
bv 


with    v=  1-t  =  a'+b't;     a' =  1,  b' = -1 
u=l  +  t  =  a  +  bt;     a  =  1,  b  =  1 
k  =  ab;  -  a'b  =  -1-1  =  -2 
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cos  vi  j^        r  .    /mfos  vl 

—  =    /l  -  tz  +  2  tan        / ,  _  J 


cos  V2 


/1+t 


cos  V2 


By  the  identity 


cos 


ix=7t 


+  cos  X 


cos  X 


=  tan  2  (n  -  x) 


--1  J^^  -  K«  -  0 


A-40: 


I16  =  sin 


v^  -  sin  v2  "  Vj^  +  V2 


For  n  =  2:    fn(t)  =  2t  /l  - 12 


I 


cos  Vl        rj-r^ 
26  =  2j  t./4-^-dt  -  2 

cos  V2 


1  +t 


cosv!  t(1.t) 


cos 


^2 


n-t< 


dt 


by  161,  Ref.   (3) 

xdx 
J  /ST? 

by  171,  Ref.  (3): 


/a2T? 


T    x2dx 


x    /— 9"     2  ,  a^ 
=  -  77  /a^-xz  +^5- 


.    -I  x 

sin  x  — 


with  x  =  t,  a  =  1 


^^{(-l^/T^-lsin-lt}-^! 

cos  V2 
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*26  =  2^s^ri  vl  "  s^n  V2)  +  cos  vl  s^n  vl  "  cos  ^2  s"1  ^2  " 

■G±viX?±v2) 


A-41: 


L^  =  2  sin  V2  -  2  sin  v^  +  2  (sin 2  v^  -  sin 2  V2)  -  (±V^)+  (±V"2) 

Note:    Sign  is  negative  for  v    in  first  or  second  quadrants  and  positive  for  third 
and  fourth  quadrants. 

For  n  =  3:    fn(t)  =  4t2  A~t2  -  A  -t2 
Following  the  method  used  in  l26: 
cos  vi 


I36-4 


cos 


t2(l-t)dt 


cos 


vi 


cos  V2    /TTT 


„cos 


vi   t3dt      .  pcos  vi  t2dt  -  rcos  vl  /rr 

=   -4  J  +41  _.  ■ 

cos  V2      A  _  t2  cos  V2     A  "  t        cos  v^    /1+t 

By  175,  Ref.  (3): 

with  a  =  1,  x  =  t,  the  first  integral  yields: 

-  -o  (sin^    vi  -  sin3    V2)  +  4  (sin  v^  -  sin  vp 
The  second  integral,  from  L^  yields 


-2 


126    -  4  (sin  vi  -  sin  vp 
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The  third  integral  =  -  I  I16J 
A-42: 

x36  =  §  {sin3  v2  "  sin3  vl}  "  2  [l26]  "  [J16] 

Considering  1^  and  In^: 

These  integrals  must  be  solved  by  use  of  Elliptic  integrals  and  Jacobian 
Elliptic  Functions. 

The  solutions  for  these  integrals  can  be  found  by  transforming  them  into 
algebraic  form  as  was  done  for  the  integrals  1^,  In3»  In5»  an<^  *n6*     Upon 
transformation  it  is  found  that  they  all  are  of  the  form  solved  in  series  235  of 
Ref.  (5).    In  utilizing  these  solutions,  the  parameters  will  have  the  following 
values: 

a  =  1,  b  =  -1,  c  =  -I,  yp  =  cos  vp 

k2  =  a^b_  m  _J_.      /2_  2  =  ill 

K         a-c        <t+l  '    K     "  x    K         -L+l 


2  2  2  / ;    , — i"      i    ,  i  '2  2^ 

g  =  =  ;  -^  =  /£+  l  ;   l+k^  =  j-r 

/a^c        /TTT     kz  ^+  1 


_!     /(a-c)(yp-b)  ^    /«,  +  l)(l  +  cos  vp] 

L"    V(a-b)(yD-c)     "  Sin      V  2(t+cos  uJ 


CPp  =  amuP  =  Sm"IV(a-b)(yp-c)     =sin"V 2(t+cos  ^ 

Up  =  F(cpp,k)  =  Incomplete  elliptic  integral  of  the  first  kind. 
E(cpp,k)  =  E(up)  =  Incomplete  elliptic  integral  of  the  second  kind. 
TT(cpp,a2,k)  =  n(up,cr)  =  Incomplete  elliptic  integral  of  the  third  kind. 
Vl   =  TT  -($+©)  J   V2  =  tt  -  (cp  -  ®) 
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Note:    Due  to  the  restriction  on  this  group  of  integrals  that 
(a  ^  y  >  b  >  c)   yp  >  b  =  -1 

Considering  Ini: 


Inl  « 


'^2    cos  nv  cos  v  dv 


vi      [£+  cos  v] 


372" 


*-J 


cos  vl  gn(t)  t  dt 


COS   V9     (-t+t)    A  "t2    A+t 


"nl 


=  r 


cos  vi 


gn(t)  t  dt 


cos  V2      (UtJ/TT/T+t  /l+t 


For  n  =  0:    gn(t)  =  1 


I 


cos  V2  r     t     -| 

COS   Vo 


dt 


/1-t  /1  +  t  /I+t 


Considering  the  bracketed  term: 


«  cos  vi 


L01 


dt 


cos  v2    /T1!  /t-(-l)  /t  -  (-£) 


cos  vi 


dt 


-if 

Jcos  V2   Ct-(-t)]  /r^t  A-(-l)  /t  -  R5 


By  235,  00: 


dt 


/(a-g)(t-b)(t-c) 


=  gF(cpp'k) 


with  a  =  1,  b  =  -1,  c  =  -I,  the  first  integral  in  IQ1  => 
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2 

—    [pfa.k)  -  F(cp2,k)]  =7y=  [u1-u2] 


By  235.01: 


r-Yr 


dt 


& 


o       (t-c)  Aa-t)(t-b)(t-c) 


^)EVk) 


with  a=l,  b  =  -l,  0=-^  the  second  integral  in  Iqi  =) 


21 


(l-l)/Z+T 


[E(cpi,k)-E(cp2,k)] 


A-43: 


101  =  ~/TTT  (U1  "  U2  "Th  [E(ul}  '  E(u2>]} 


JT+T 
For  n  =  1;    gn(t)  =  1 


Lll 


pCOS  vi 

cos  V2 


ff+ty 


dt 


/FT/T+t"  /I+t 


Considering  the  bracketed  term: 

t2         .  [(£+t)-l,]2     .  (£+t)2  -2t(l,+  t)  +  l,2 
(I+tJ  (£+t5  ft+tj 


a+t) 


=  <l  +  t  -  2£  + 


tf 


-t  +  t 


=   t-<t+- 


^ 


1+  t 


t-^fl 


£T+T> 


d 


comparing  with  Iq,: 


cos 


[H  =  J. 


vl 


t  dt 


-I    I 


cos  V2    /T^ryr1  pH  /t-(-'t) 


01 J 
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By  235.  16: 


I 


yp 


tm  dt 


b        /(a-t)(t-b)(t-c) 


=  gbmJ 


C^      ck2       2    N 
^     ^l-„sn2u^ 


(1  -  k2  sn2u) 


du 


m 


with  m  =  1,  a  =  1,  b  =  -1,  c  =  -I 


gbm=li) 


t+t    /ttt 


The  integral  in  I-q  becomes 


PU1  (1  -  £k2  sn2  u)    du    +2       ru2    (l-tk2sn2u)    du 
o      (1  -k2  sn2  u)  JUT  Jo       (1  -k2  sn2  u) 


Jwi 

By  340.  01  with  a2  =  k2,  a^  =  tk2 


rilipgjLdn-4-   r(k2-^2)n(cp,k2,k)  +  tk2u] 
J  1  -  k2  sn2u  k^      L  J 


=  [(1  -«t)n(cp,  k2,  k)  +  <tu] 


A -44: 


111  ■  ^=*   fndi.k^k)  -  n($2.  k2,  k)l    -=L(Ul  -u2)  -  ifl^] 

Where,  by  111.  06,   Ref.   (5): 

E(cp,  k)  -  (k2  sin  cp  coscp//l  -  k2  sin2  cpj 


n(cp,  k2,  k)  = 


For  n  =  2:    gn  =  2t2  -  1 


k' 


cos 


I 


^    r2t3-ti 


dt 


=  r 

21       Jcosv2      L<^J     /(i-t)[t-(-l)][t-(-t)] 
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Considering  the  bracketed  term: 


2t3-t   =  2[q+t)  -  if  -  1(1 +  t)  -  f] 

C-t+t]  i  +  t 


2(£-K)3  -  6(£+t)2t+  6(£+tK2  -  2t?  -  jl+X.)  + 1 


=  2(t+t)2  -  6^+t)  +  6t2  -  1  -  \jfc+[f  ] 


2t2  +  4£t  +  2t2  -  6£2  -  6U  +  6£2  -  1  -  [ffi"i  j 


=  2l2  -  2£t  +  2t2  -  1 


\-2l3-l  1 

L(^  +  t>  J 


=  2t2  -  21  (  t  - 1  + 


t"t+CLTtT^  +  7I+t 


Comparing  with  1^ ^  and  Iq^: 


I21  =  2 


I 


t2  dt 


/(l-t)[t-(-l)][t-(-t)] 


-  2*[ln]  -  [ioi] 


Again  citing  235.  16,  with  m  =  2 


I 


t2  dt 


A -t)[t-<-i)][t -<-*)] 


-i 


Ul  9         9      2 

(1  -  UT  sn2ur  du 


(1  -k   sn    u) 


■"I1 


U  9       9       2 

(1  -  £kz  snz  u)   du 


(l-k2sn2  u)2 


By  340.  02:    with  or  = 


u2        2 

k  ,  ax 


=  -tk" 
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f  a-^^Vdu         1  f  A4U  +  2^,2^2  .  ^2)v    +  (k2  .  ^2)2       1 
J   (1-k2  sn2u)2  k41"  2J 

By  315:    (since  a2  =  k2,  Vj  =  I2j) 

Vx(u)  =  ~j  [E(u)  -  k2  sn  u  cd  u] 
k' 


V9(u)  =  -i-r  [2(2-k2)E(u)  -k/2u]-k2snucdu(k/2nd2u+4-2k2)] 
3k' 4 


4(2 -k2)  V0-3V,  -k2  snu  cdund4u 
V3(u)  =  1 L— g 


5k' 

Note:    Vp  V2,  Vo  are  used  to  evaluate  both  I2p  lop  and  I04.    They  are  most 
readily  applied  by  evaluating  them  and  substituting  in  their  values. 

J  (Wk2Sn2U)2       =  W*  +  2^  ~  *)  V2(u)  +  ( 1  - D2  V2(u)J 
(1  -  £k    sn   u) 

Combining 
A-45: 

I21  =  2*,2<u1  -  u2)+  441-*,)  [v2(ui)  -  Vx(u2)] 

+  2(1  -  I)2  [v2(Ul)  -  V2(u2)]  -  21  [ln]  -  [lQ1] 

For  n  =  3:    gn(t)  =  4t3  -  3t 

i      =   ,cos  Vl  -  4  4  .  3t2  -J g 

31       4os  V2    ^     {l+t)   J  /(l-t)[t-(-l)][t-K)] 
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Considcring  the  bracketed  term: 


4t4  -  3t2     .  t2(4t2  -  3) 

(i+t) '        IX+tj 


t2  =  1(1  +  t)  -  l\2  =  (1  +  t)2  -  2U,l  +  t)  +  I2 


2,t.2 


t  (4t"  -  3)     _   2  rA(J  ,  tx     o .  ,  4^-31 


=  4t3  -  4£t2  +  (4^2  -  3)  U+  t)  -  2l\  +  4 


w 


4  2 

4t3  -  4£t2  +  4£2t  -  3t  -  4l3  +  31  +  4J 


T^+ty 


2-e?- 


-  4t'  -  2i  {2^  -  2«  +  2t2  -  1  -  [^J}  +  4-Ut  ^y 


3.,  -i 


-i.2 


=  4t3-2t{2t2-2«+2t2-l-2^i]}  +  3(i-t^)-2t[l-1;r 


comparing  with  I^2»  Iji*  and  I 


10 


COS  V 


131  =  4  L  *  a:  mMkm  -^  ■ M  &  J  -  ftJ  - 2* 


oi  J 


considering  the  first  integral: 


By  235. 16,  with  m  =  3: 


rUl  (l-£k2  sn2  u)3du     _4  TU2    (l-£k2sn2u)3du 
°o       (l-k2sn2u)3  Jo         (l-k2sn2u)3 
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By  340.  03,  with  a2  =  k2,   ctj2  =  <tk2: 


p  (l-^2sn2u)3    du  m    1     1-^6     +  3t2k4  (k^Jc2)  y    ,  g  ^2(k2  .  ^2 
J      (l-k2sn2u)3  k6    L  *  2 

+  (k2  -  ^2)  v3] 

J(l-^c    sn    u)3    du  =  r/t3u+3/t2(1_/t)Vi+3/t(1_/t)2V2  +  (1_^)3  y    1 
(1  -  k2  sn2  u) 


A-46: 


I31  =  4t3  (ux  -  u2)  +  12 12  (1  -  t)  [ViCuj)  -  V^) 


+  12^(l-<t)2[v2(u1)-V2(u2)] 

+  4  (1  -  if  [v3(Ul)  -  V3(u2)]  -  ^[l^]-  3[lu]  -  2^[l01] 


Considering  1^: 


\>r 


I 


sin  nv  sin  v  dv 


n4      j         P ,  _  r .  -,1 

v^     LI  +  cos  v]  11+  cos  VJ2 


cos  vi        fn(t)  dt 


Jcos  V2      [1+t]  [^+t]2 

Series  235,  Ref.  (5)  will  again  be  used  to  solve  these  integrals. 

For  n  =  0:    fn(t)  =  0 

A-47:     I04  =  0 

For  n=  1:      fn(t)  =  ^(1 -t2 


lU  = 


dt 


.cos  vip/f^T- 

=osv2  L    (1+t)J"[^]J' 
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considering  the  bracketed  term: 


A  -t2       /!Tt  /l  +  t         / 


1-t 


l+t  l+t 


/l  +  t 


L14 


•cos  vi         /FT  dt 

cos  V2    y(t-(-l)][t -(-£)] 


By  235.  07: 


fV^-b^V-c)   dt  =  (a-b)g|oUPcd2udu 


with  a  =  1,  b  =  -1,  c  =  -I, 


L14 


Ui 


cdz  u  du 


ru2 


Sl+1   Jo  Jl+1    o 


cd2  u  du 


By  320.  02: 


I 


cd2  u  du 


—c)  [u  -  E(u)  +  k2  sn  u  cd  u] 
kz 


A-48: 


I14  =  2/T+T    ui  -U2  +  E(u2>  -  E(u1)  +  yxj(snuicdu1  -  sn^cd^)! 


For  n  =  2:    fn(t)  =  2t  7(1  - 12 

cos  vi  i - 

r2t/l-t2l 


L24 


dt 


cos  v2 


L-d+tTJ^T^r 
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Considering  the  bracketed  term: 


2t /l  - 12       ItJITt        nn.^/T^t       0  /l  -t 

— 7V+ 1\ : 2(l+t)    ■    •■■-  -  2    - 

(1+t)      /r+r  /r+t      /ttf 


=  2 /T+t  /T-t  -  2 


1-t 

/r+r 


Comparing  with  I,  A 


hi-2! 


cos  vi 


COS 


^2 


T+t /FT 
/I+T 


dt  -  2 


M 


By  235.  09: 

ryp 
4> 


(a-t)(t-b)    _ 


"IF^T 


=  (a-b)(b-c)gk' 


sd2  u  cd2  u  du 


With  a  =  1,  b  =  -1,  c  =  -l,  the  integral  in  I24  => 

Ml 


2t2(l~ x)  (  T^Xttt)] t  f1  sd2  u  cd2  u  du "  [  ' sd2  u  cd2  u  du] 


By  361.27: 


[sd2  u  cd2  u  du  =        .1    .    |"-2k'2u  +  (l  +  k/2)  E(u) 
J  3k4k/2   L 


+  k2  (k/2  nd2  u-  1  -k/2)  sn  u  cd  uj 


since: 


1      .  (£+l)2        l  +  k/2         21         k2  2 


k'' 


t-i  '  k,2    <wy 
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Nd2  ucd2udu  =  ^iy 


12 


-2u  +  -^-   E(u)  + 


l-l 


£+1 


(nd2  u  -  .  _  v  J  sn  u  cd  u 


A-49: 


Z24  ■  3  (* -!)/£+!  {2u2  -2u1+-^-[E(u1)-E(u2)] 


2    r,  j2 


+z+r  L(nd  U1  "i^rjsnuicdui 


'Cnd2u2  ■qFr)snu2cdu2]}-2[Ii4] 

For  n  =  3:    ^(t)  =  (4t2  -  1)  /l-t2 


L34 


COS   V 


COS 


Y2 


1  r^t2-!)/!-*2] dt 


L  (T+i 


(t  +  t)2 


Considering  the  bracketed  term; 


(4t2  -  1)  /l-t2     ..  -    (2t  A+t2)       /l-t2 


TT+^ 


=  2t 


(i+t)       "TT+tT 


Substituting  2t  =  2(1  +t)  -  2  => 


4t 


/T72*   -*£ 


/TT2        /T"? 


(i+t)     "Tr+ty 


Comparing  with  I41  and  I42: 


cos 


cos  V2       /£+  t 


[24]  "  \ 


T 

X14J 


with  t  =  (t  +  t)  -  t,  the  integral  in  I43  becomes 
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COS   Vi  COS   V)         /       i 

4  J  Jd-t)[t-(-i)][t -(-*)]  ft  -  U  J        »  /<1  -OCt^-DJ  dt 

cos  V2  cos  V2  /t  -  (-<l) 

The  second  integral  is  the  same  one  as  was  solved  in  I42 

_4l  rcos  Vl  Ai-t)Lt-(-i)]'  =  _ur  h*  +  2Ii4 1 

Jcos  V2  A  -  R3  2 

Considering  the  first  integral: 
By  235.  14: 

pV  u 

Jo       ,y(a-t)(t-b)(t-c)dt  =  g(a-b)(b-c)2k2  J       sd2ucd2  und2  udu 

with:    a  =  1,  b  =  -1,  c=  -I,  the  first  integral  =) 

4^ ~^)(2>('C'"1)2CtTT-)[  sd2ucd2und2udu  -  J       sd2ucd2und2udu] 

By  361.  18: 

|sd2ucd2und2udu  =  -~  [-I2  +  (l  +  k/2)  I4  -  k'2  Ig] 

where     I2  =  V]_ 

I4  =  V2  as  defined  in  I21 

I6  =  V3 
A-50: 

I34  =  8(4  -  l)2/I+T{Vl(u2)  -  Vj  (ux)  -f-2^-  [v2(ui)  "  V2(u2)]  +  (±£) 


[V3(u2)  -  V3(Ul)]  }  -  2(4+1)  [l24]  -  (44+1)  [l14] 
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APPENDIX  C 

COMPARISON  OF  RESULTS  WITH  KNOWN  RESULT  FOR  A 
COMPLETE  SHROUD 


By  the  Kutta  Joukowski  Law: 

(Equation  P-2)   Fr(9)  =  -p  V^  (1,9,0)  r(Q) 

Taking  the  total  lift  in  the  9  =  0  direction: 

TT 

T(0)  =  -P  J      V^a,  9,  0)  r(Q)  cos  9  ri/4  d9 


-TT 


In  the  absence  of  a  propeller 


Vxp  =  Vw 


;  T(0)  =  -pV^  r1/4  J      r(e)  cos  9  d9 


-TT 


By  A-23: 

CO 

r(Q)  =    ^  Ancos  (n9) 


n=o 

CO  JJ 


T(0)  =  "PV^i/a    /  cos  n9  cos6  d6 

n=o    -n 


Because  of  orthogonality. 


.TT 

cos  n9  cos  9  d9  =  0;n  ^  1 

-IT 


.'.  T(0)  =  -pVror1/4  Ax  J       cosz9  d9 

-TT 
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A-51: 

T(0)  =  -pV^r^Ajir 

To  find  the  side  thrust  on  a  complete  ring,  in  the  absence  of  a  propeller, 
we  need  to  evaluate  Al. 

For  a  complete  shroud,  of  constant  chord,  at  an  angle  of  attack, the  side 
thrust  is  independent  of  the  section  shape  (Ref.  6).  We  will  consider  the  case 
of  a  complete  shroud  consisting  of  a  frustrum  of  a  cone  as  shown  in  Fig.  A4. 


*x 


FIGURE  A4 


The  shroud  has  a  chord  of  L,  a  radius  at  the  1/4  chord  point  of  r1/4,  a  con- 
vergence angle  3  and  is  set  at  an  angle  of  attack,  in  the  9  =  0  plane  of  ct(0)  =  ccQ. 

C'  ($,X)  =  0 

a($)  =  -$  +  Oq  cos  $ 


•71- 


By  Pi: 
A-52: 


<lr< 


By  A-26: 


=  rV   [+p  -  a0  cos  §] 


(l,cp,\3/4) 


CD  6 

5rs  =  ?    I   I       Cnm  *a 
n=o  m=l 


m 


In  evaluating  the  integrals  Inm  we  must  abide  by  the  restriction  that 

Vp    \>2    f   TT 

where 

VI    =   TT  -  ($  +  ©) 
V2    =   TT  -  (§   ~  ®) 

Since  ©  -  TT,  for  a  complete  ring, 

$  7^  TT 

V2  =  2tt  -  $  ;  Vo  "  vi  =  2  tt 

Note:    In  evaluating  the  elliptic  functions   E(cpp,  k)  and  F(cpp,k),  the  following 
relations  are  used: 

E(-cp,k)  =  -E(cp,k) 

F(-cp,k)  =  -F(cp,k) 

E(mTT  ±cp,k)  =  2m  E(k)  ±  E(cp,k) 

F(mn  ±  cp,k)  =  2m  K(k)  ±  F(cp,k) 
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where: 

E(k)  =  Complete  elliptic  integral  of  the  second  kind. 

F(k)  =  Complete  elliptic  integral  of  the  first  kind. 
Determination  of  cpi  and  cp2  from  v^  and  \k: 


/  (£  +  1)  (1+COS  Vp) 

^P  =  sin"    V        2«,+cos  vp) 


2  sin2  $p  (-t+cos  vp)  =  (-1+1)  (1  +  cos  vp) 

=   i,  +  t  COS  Vp  +  1  +  cos  Vp 

(£+cos  vp)  (2  sin    cp«  -  1)  =  1  +  K,  cos  vp 

1  +  I  cos  vp 
■cos(2(^P)    =    t  +  COSVp 


when  vp  changes  by  2tt,   |2cpp|   changes  by  2tt.    To  find  the  sign,  assume 
1=2. 


(1)  let  vp    =  0 


-cos  2cpp     -  2  +  1  =  1  ;    cos  2cpp     =  -1   ;  cpp     =  -j- 


(2)  let  vp     =  tt/3 


-  cos  2cpp  =  =  0.  8 


2+1 


2 


2cpp2  =  tt  -  37°  cpp2  cpp2  =  5  "  18° 


(3)  let  vp3  =  rr/2 

-cos2cpp  =  2  5  cos  2cpp  =  "2 
2cpp  =  2  tt/3    cpp  =  n/3 
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(4)  let  vP    =  n 


-cos  2cpp4  =  -l/i  =  -i 
cos  2cpp4  =1    cpp4  =  0 

when 

vp  changes  by  +2  Tr 

cpp  changes  by  -tt  as  shown  below: 


-  la. 


E(cpP  -  n,kf-  E(cpp,k)  =  -2E(k)+  E(cpP,k)  -  E(cpp,k) 

E($p  -  n,k)  -  E(cpP,k)  =  -2E(k) 
Similarly 

F(cpP  -  TT,k)  -  F(cpP,k)  =  -2K(k) 

Since  in  this  case 
v?  -  V]_  =  2rr 
$2  -  cpj  =  -TT       $2  =    $j_  -  n 

E(u2)  -  E(Ul)  =  E(9l  -  n,k)  -  E(cpi,k) 

E(u2)  -  E(Ul)  =   -2E(k) 
Similarly 

u2  -  Ul  =  -2K(k) 


•74- 


By  A-43: 


L01 


/T+ 


=  {+2K(k)-T^rE(k)} 


where 


I  = 


\2  +  2 


'    l~l  ~  T 


,2  2 

k     = 


£+1 


=     \z  +  2  +  1 


(X2  +  4) 


k'z  =  1-k*  = 


2  _  t-1 


1+  1 


A? 


(t+i)  =  ii+i.j/E+r=/^±± 


-t,  X2  +  2 


T=Y       ^J 


i/T    r  r\2  +  2 


01    /7T7  l  v   ,2 


[-(}-¥- )E(k)  +  K(k)] 


I01  =  4/2  [-  ( 


X2  +  2 


X2S? 


+  4 


'  A2  +  4 


A-53: 


L    Com  W  "  C01  ^1 
m=l 


_      "Ao     V_(       \2  +  2    n  \K(k)     1 

"1/4     L    ^2X/?71J     {)      2/Z711 


-75- 


For  n  =  1: 


X12  "  X13  "  X15  "  ° 


By  A-44: 

T       -  2(t-l)    r+2E(k)  1 


Since 


[+2K(k)]-t[Ioi] 


.2 


En  =  — = —  (+2(^+1)  E(k)  +  (37^)  E(k)  -  2t  K(k)  -  2£K(k)} 


ln  =  ^r  (+[G^r)  +  <^  +  *>]  E<k>  "  2^K^>} 

By  A-48: 

I14  =  2  /T+T  [+2  K(k)  -  2E(k)] 


hA=-4=z    [-(*+!)  E(k)  +  a+l)K(k)] 


L14 


By  A-40: 


I16  =  +  2tt 
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6 

I 

m=l 


L    ClmIlm  *  C11I11  +  C14I14  +  C16I16 


.  XAicos^  pT   ,T     /2  T  -1 


l^-7m 


8/2nr1/4 


4  r      y    *2 


|+(^r)E(k)-(t-i)K(k)} 


.  .,  X2+2    .',     .     .X2     .    .  ,  ,     .  X2  +  4  .    £2      ,  Q,2  +  2)2 


In  +  I       =  -ML.  {+{X    +22)       E(k)  -  \2  K(k)} 

A2+4  c       r  J 


A-54: 

6 


_,  nrl/4       L     4XA2+4  4  A2 +4         *J 


Similarly,  it  can  be  seen  from  the  forms  of  Cnm  and  Inm  that  for  a  com- 
plete shroud: 

£  A    cos  2  $ 

I    <W2m--      nrlM      [H^k)] 
m=i 

and  in  general: 

A-55: 

6  (-l)n+l  a    cos  n$ 

>      CL^L,™  = [Hn(X,k)] 

l_j       nm  nm  tt  r,  /4  n 

m=l 
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where  Hn(X,  k)  is  some  function  of  \  and  k  to  be  evaluated.    Substituting  A-53, 
A-54,  and  A-55  into  A-26,  and  equating  A-52  to  A-26: 


TTrl/4  L    ^  2\f\2  +  4'/  2  A 


XK(k)    1 


+  4 


A.  cos 


2  +  2)2    "> 


nrl/4      L     L4\A2  +  4J  4  A2 +  4       4J 


n=2 


"     (-l)nAncosn$ 


n  r 


1/4 


[Hna,k)] 


Equating  coefficients: 


A     = 
o 


-V-nrl/4P 


r  xz  +  2    n 


LA 


2\  A2+4 


y 


E(k)- 


XK(k) 

2  A2 +4 


Al  = 


-V_  tt  r,  ,„  a 

00        1/4    o 


1     L4\A2T4    J  4^274        4    ; 


An  =  0,  n^2 


Substituting  Aj  into  A -51: 


T(0)  = 


+PV^r2/4aQ 


{  +  [J^]E(k)-^^+£} 
L     L4\A2+4J  4A2  +  4        *J 
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Recalling: 


rl/4 


\2  +  4 


■; 


and  utilizing  a  lift  coefficient  analogous  to  that  used  in  Ref.   (2): 


CL  = 


T(0) 


T(0) 


(p/2)v2Lr1/4         pVl\r21/4 


+ 


r? 


cu 


{♦[ 


(X2  +  2) 


2    -, 


E(k)  -    A<k) 


4X A2  +  4  J  4  A2  +  4 


3-} 


A-56: 


+  4tt- 


^[t^]™-2*^"} 


/? 


+  4 


A2 +  4 


This  result  is  equivalent  to  that  obtained  by  Weissinger  in  Ref.   (1)  and,  as 
shown  in  Figure  A -6  is  very  close  to  the  result  obtained  from  the  more  exact 
theory  by  both  Weissinger  and  Morgan  in  Ref.  (2). 
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